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ABSTRACT 

Precise numerical standard solar models (SSMs) now agree with one another and with helioseismo- 
logical observations in the convective and outer radiative zones. Nevertheless these models obscure 
how luminosity, neutrino production and (/-mode core helioseismology depend on such inputs as 
opacity and nuclear cross sections. Although the Sun is not homologous, its inner core by itself is 
chemically evolved and almost homologous, because of its compactness, radiative energy transport, 
and ppJ-dominated luminosity production. We apply luminosity-fixed homology transformations 
to the core to estimate theoretical uncertainties in the SSM and to obtain a broad class of non- 
SSMs, parametrized by central temperature and density and purely radiative energy transport in 
the core. 
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1. Homology and the Solar Core 

After more than three decades of nuclear cross section measurements, opacity calcula- 
tions, and detailed computer evolutionary calculations, standard solar models fSSMs) with 
the same inputs now agree in their neutrino flux predictions to within about 1%. The the- 
oretical models are now also consistent with precise p-mode helioseismological observations 
of the Sun's outer radiative zone x = t/Rq — 0.26 — 0.71 and convective zone x > 0.71. If 
the ^-mode helioseismological oscillations have sufficient amplitude, their observation can 
be expected soon to calibrate the solar inner core, where the thermonuclear luminosity 
and neutrino production take place. While necessary in the complex convective zone and 
justified by the precise helioseismological observations, the complexity and numerical form 
of precise SSMs obscure the simplicity of the solar core and the determinants of solar neu- 
trino fluxes. In order to understand standard and non-standard solar models, we return 
to the homology methods of Schwarzschild (1958) and Iben (1969 and 1991) used before 
the advent of fast computers, but with three new features. 

Castellani et al. (1993) have found that changing input parameters by factors as 
large as two leads to only homologous changes over 60% by mass of the Sun. In this 
paper, we explain this remarkable homology and demonstrate that, while the entire Sun is 
certainly not homologous, the core is homologous enough to be parametrized by its central 
temperature, Tc and density pc- This (Tc, pc)-pai'ametrization subsumes all astrophysical 
effects of opacity, composition, and the ppl nuclear cross section factors S\\ and 
(the pp and ^He — ^He reactions) into the two parameters (Tc, pc), one representation of 
the central boundary conditions of solar structure. Indeed, any standard or nonstandard 
solar model that depends principally on radiative energy transport can be parametrized 
by (Tc, Pc) and the remaining nuclear cross section factors 534 and Sn (the ^He — ^He 
and p — ^Be reactions) (DegliTnnocenti 1994). Even solar models with a non-standard 
low opacity or low metallicity Z are all essentially parametrized by {Tc, pc) or by Sn, 
the principal cross section factor determining Tc. (See Figure 2 of Hata et al. (1994) or 
Figure 2 of Hata and Langacker (1995); see also Hata (1994).) As the pc dependence of 
the neutrino fluxes is weak (see section 4 below), Hata and Langacker (1994) were able to 
show that the 0.7% theoretical uncertainty in Tc, together with the remaining nuclear cross 
section uncertainties, provide the same theoretical neutrino flux and rate uncertainties and 
correlations that Bahcall and Ulrich obtained from 1000 Monte Carlo SSM simulations. 
(See Figures 2-4 and 6-8 of Hata and Langacker (1995).) 

The (Tc, Pc) parametrization allows analytic estimation of the logarithmic derivatives 
Pj{i) = d\n(j){i)/d\nSj of the principal neutrino fluxes (f){i) with respect to input parame- 
ters Sj, which Bahcall and Ulrich (1988) obtained from 1000 Monte Carlo SSMs calculated 
with sma// changes in input parameters. Because homology makes these logarithmic deriva- 
tives constants, from any precise SSM, we can not only estimate theoretical uncertainties, 
as did Bahcall and Ulrich, but we can now also extrapolate to non-SSMs, so long as the 
energy transport is primarily radiative. 

Our application of homology differs from earlier ones in three ways: (1) We apply 
homology only to the solar inner core, not to the whole Sun; (2) we do not assume p ~ T*^ 
or any polytropic relation; and (3) we use homology to derive the dependence of core 
temperature and density on opacity, nuclear energy generation, and mean molecular weight, 
at fixed luminosity, instead of the dependence of effective temperature on luminosity, with 
fixed opacity and nuclear energy generation (Cox and Giuli 1968). 

After accounting for the different energies released when pp, ^Be, ^B, and CNO neutri- 
nos are produced, the known solar luminosity Lq fixes the total photon energy production 
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and constrains the neutrino fluxes through the nuclear reactions: 

(f){pp) + (0.967) (piBe) + (0.743) (p{B) + (0.946) (j){CNO) = 6.48 x 10^° cm-^g-l , (1) 

for the four principal neutrino fluxes. (See Appendix; Castellani et al. (1993); Bludman 
et al. (1993); Hata and Langacker (1994).) For the SSM, the three pp branches /, //, /// 
terminate in the ratio 83.7% : 16.3% : 0.02%. Differentiating this constraint and assuming 
these termination ratios continue to hold approximately, we have the constraint 

f^iipp) + (0.079)/?i(Se) + (0.000071)/?i(S) + (0.0145)/?i(C7VO) = (2) 

on the logarithmic derivatives of the principal neutrino fluxes with respect to any input 
parameter Sj. The logarithmic derivatives obtained by Bahcall and Ulrich from their 1000 

Monte Carlo SSMs satisfy this sum rule (Bahcall and Ulrich 1988; Turck-Chieze, 1988; 
Bludman et al. 1993; Hata 1994; Hata et al. 1994; Hata and Langacker 1994, 1995). 
For models that depart greatly from the SSM, the three termination ratios in general 
change significantly from the values used here; in that case, the new ratios can be used as 
the new constraint, for small deviations from the new SM. The numerical coefficients of 
equations (1) and (2) also change. 

The central temperature and density are outputs characterizing solar models. Along 
with the nuclear cross sections and chemical composition, given by the vector X of element 
abundances by mass, they determine the neutrino fluxes. The logarithmic derivatives of 
these fluxes with respect to central temperature, a{i) = d\n(f){i)/d\nTc, satisfy 

a{pp) + (0.079)a(Se) + (0.000071)a(S) + (0.0145)a(CiVO) = , (3) 

and are functionals of temperature, density, and composition, that can be approximated 
by power laws in temperature, if the weak pc dependence is ignored. If aijip) — —1.2, then 
a{Be) = 8±2, a{CNO) = 34±9, is consistent with this constraint. Once the ppll, ppIII, 
and CNO cycles are included, the luminosity constraint prevents the solar core from being 
strictly homologous and induces a small pc dependence. 

2. The Present Sun 

The luminosity- and neutrino-generating core and the outer radiative and convective 
zones are almost decoupled dynamically. For this reason, the solar model inputs of luminos- 
ity Lq and radius Rq almost separately determine the two outputs, the helium abundance 
Y and the convective zone mixing length (Bahcall and Ulrich 1988; Turck-Chieze et al. 
1988). This permits us to ignore the convective zone, for which accurate opacities and 
detailed numerical models are needed, and to concentrate on the outer radiative zone and 
the inner core. The radius Rq being irrelevant to the inner core leaves only the luminosity 
Lq and the mass Mq as fundamental parameters. 

The equation of hydrostatic equilibrium is 

-dP/pdr = g = Gm/r^ (4a) 

or 

= l/Xp = g/{P/p) , {4b) 

where Ap = {—dliiP/dr)~^ is the pressure scale height. We deflne the stiffness (effective 
polytropic exponent) T = d\nP/d\np = 1 + (l/?^eff) effective polytropic index n^Q = 
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din p/dln{P/ p) and write the equation of state as P/p = ^{T/p){l + D), where 3? is the 
gas constant, p is the mean molecular weight and D includes all corrections to the ideal 
gas equation of state. Then 1 — F"^ = dln^P/ p)/dliiP = V — V^, where the thermal 
gradient V = dlnT/dlnP and chemical gradient = (iln(///(l + D))/dlnP. If n^Q and 
p/{l +D) were constant, the Sun would be a polytrope of index n^Q and thermal gradient 
V = l/(ngff + 1). These conditions obtain in the convective and outer radiative zones 
of the Sun, but not in the chemically evolved and inhomogeneous inner core. The core 
structure is usually obtained by evolutionary models that depend on the initial relative 
metallicity Z/X for the proto-Sun and the present age. 

Figures 1 and 2 show the gradients V, 1 — F~^ and F as function of the dimensionless 
radius x = t/Rq and included mass m/MQ, derived from the SSM of Dearborn 1994. 
(This model provides dense enough output to allow us to plot the P, p, T logarithmic 
derivatives and agrees well with other SSMs, such as the SSM of Bahcall and Pinsonneault 
(1992) with helium diffusion.) The Sun's convective zone is an n^Q = 3/2 polytrope; 
the outer radiative zone, approximately an n^Q = 4.3 polytrope. The core, chemically 
evolved and inhomogeneous, contains the Sun's luminosity production and the majority 
of its mass. This complex structure prevents the entire Sun from being homologous, even 
though luminosity production and opacity are approximately power laws in individual 
zones. 

The simplicity of the passive outer radiative zone governs any matter-amplified neu- 
trino oscillations which may take place there. The concentration of mass and luminosity 
production interior to this zone makes it approximately an UqQ — 4.3 polytrope, with 
density scale height Xp = TXp very nearly constant at Rq/ 10. 5. With this scale height, 

if neutrinos of energy E and mass squared difference Am? oscillate with vacuum mixing 
angle sin^ 6, the adiabaticity is 

A = 5.3xl0'^{Am'^{eV^)/E{MeV)){sin^2e/cos2e) , (5) 
and the jump probability Pj = exp{—E]\iji^/E), with 

EjsiA = (71^/2) AE = 3 X 10^ Am2(sin2 2^/cos2^) . (6) 

The constant density scale height of the n = 4.3 polytropic outer radiative zone is the 
principal property of the Sun effecting the Mikheyev-Smirnov-Wolfenstein (MSW) neutrino 
oscillations. The best small-angle MSW solution to the solar neutrino observations gives 
^NA ~ 13 MeV, so that neutrinos of energy E=10 MeV oscillate non-adiabatically 
with A = 0.3 and pp neutrinos of energy E=0.3 MeV oscillate adiabatically with A = 10. 

The solar inner core is reasonably inferred to have been initially convective and there- 
fore chemically homogeneous. The thermal gradient was thus initially adiabatic, V = Vad, 
decreasing since. Meanwhile, the composition gradient has been increasing from zero. 
At the present epoch, these two evolutionary changes nearly compensate, V ~ V^, mak- 
ing F pa 1, so that T//x pa constant (Figure 2) and the core is quite condensed. Together 
with the thermal gradient V fa fa 1/3, this makes T^/p constant to within 7% for 
X < 0.3, ui/Mq < 0.613 (Bahcall and Ulrich 1988). This accident of the present epoch 
makes {p/{l + D))P ~ p4/3^ 

which, in the chemically inhomogeneous core, is not an 
ngff = 3 polytrope. Since P j p T j p, constant in the core, the core structure is nearly 
an "isothermal" n^s^ — oo polytrope with T/p, rather than T nearly constant. We note for 
completeness that a long-lived convective core (such as a present-epoch fully mixed core) 
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would not follow the simple homology of this paper, as the simple chemical stratification 

T/|U ~ constant would be violated. 

Although far from polytropic, the solar inner core is almost homologous, because over 
the narrow range of density and temperature in the compact core, the nuclear energy 
generation and Rosseland mean opacity are approximated by the power laws, 

e = eoP^T\ K = Kop''T-' , (7) 

where Eq, Kq, p depend on the chemical composition. Outside of the core ranges of density 
and temperature, these homologous forms would still be valid, but with different exponents. 
(Alternative forms of energy transport, such as WIMPs, would destroy the homology 
altogether.) The core, which we define by radius x < 0.26 so as to contain 99% of the 
energy generation and 50% of the mass, has T = (7.7-15.7) x 10^ K, p = (19-154) g cm^^, 
and a small pressure P = Pc/l^ at its edge. In this range, wc adopt A = 1, u ^ 4.24 for 
the luminosity generation (dominantly but not exclusively ppl cycle) and fit the inner core 
OPAL opacity (Iglesias and Rogers 1991; Rogers 1992) by n ~ 0.43, s ~ 2.47. Our inner 
core OPAL opacity fit is close to the values n = 0.5, s = 2.5 illustrated by Figure 10.4 
of Bohm-Vitense (1992). For comparison, we also consider the Kramers opacity (Cox and 
Giuh 1968). 

3. Homology Applied to the Solar Core Only 

Homology is usually applied to zero-age main sequence (ZAMS), chemically homo- 
geneous stars to derive power-law L — M, R — M, L — T^g- relations among luminosity, 
radius and effective temperature for different stars with the same power-law opacity and 
luminosity generation. As shown in Figure 3, taken from Kippenhahn and Weigert (1990), 
these homology relations are obeyed by models and actual stars on the ZAMS upper and 
lower main sequence, which respectively have convective cores or envelopes. For these 
stars, L ~ M^-^^, R ~ M^-^^ and L ~ M^-^, R ~ M^-^, respectively. 

The Sun, however, is not chemically homogeneous, is mostly radiative, and is transi- 
tional between pp and CNO burning. Even on the ZAMS, the Sun was not homologous. 
Figure 3 shows that, for (0.3 — 1.3)M0 stars near the Sun, the exponents change rapidly 

with mass, with L ~ M^-^, R ~ M^ '^ at the solar mass. Therefore, we cannot apply 
homology to the entire present or ZAMS Sun. 

On the other hand, because the core of the Sun has approximately power-law opacity 
and luminosity generation and has a low pressure boundary condition, P k, P^/12, we can 
apply homology to the isolated solar core. We do so in a different way: instead of deriving 
luminosity, mass and radius relations for a family of stars having the same opacity and 
luminosity generation, we derive the dependence of central temperature rc(K, £, p) and 
density pc{i^i p) on overall opacity, energy generation, and mean molecular weight, for 
one star of fixed luminosity. This enables us to scale from any given SSM to models of the 
same luminosity with different input parameters. 

Assuming an ideal gas equation of state, the mass conservation and hydrostatic equi- 
librium equations scale as 

p ~ m/r'^ , P ~ /r^ , (8) 
where m the mass included inside radius r, so that 

m/r ~ P/p ~ T/p , P^/P ~ T^/P ~ {p^mf , T^/p ~ p^rr? , (9) 
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where P and are the total and radiation pressures. The first and second relations give 
the virial theorem m/r ~ T///. Any given correction D to the ideal gas equation of state 
can be included in the homology formulas by replacing with ///(I + This substitution 
ignores any p, P, or T dependence in D. Typical corrections include (Bahcall and Ulrich 
1988; Bahcall and Pinsonneault 1992): the Debye-Hiickel screening effect, contributing D 
~ —0.014; photon pressure, contributing D ~ +0.001; quantum degeneracy, contributing 
a negligibly small effect; and a hypothetical core magnetic field, contributing D ~ +0.002 
for a field strength of 10^ gauss and scaling as the square of the field strength, if the reverse 
influence of the thermomechanical structure upon the magnetic fleld is ignored. 

The equations for radiative energy transport and thermal steady state are, using equa- 
tion (4b), 

K(£/m) = ATTcG{dP^/dP) , e = M/dm , (10) 

where the Rosseland mean opacity, luminosity, specific thermal energy genera- 

tion at radius r, so that 

Ir^l/m^/K . (11) 
From equations (9) and (10), we deduce how £ scales with fx, m, r : 

SK ~ T^/P ~ {^?mf . (12) 
The quantity r^+^^s+^n ^ ^^^^^v-s-A^v+\-s+n-2 -^e eliminated to obtain 



where 



s — 3n 

a = 



v + 3X — s + 3n 
P = l + a 

(4 + 3A + 3n)(s-3n) (14a) 



7 = 4 + 3n + 
5 = 3 + 2n + 



z/ + 3A — s + 3?i 
(2 + 2A + 2n)(s-3n) 



e 



+ 3A — s + 3n 
3 + 2n 



2 + 2A + 2n 
C=l-e 

(3 + 2n)(4 + 3A + 3n) ^ (146) 
2 + 2A + 2?i 
^^(3 + 2n)(. + 3A-. + 3») 3^ 
2 + 2A + 2n 

This homology rests on equating the luminosity produced with the luminosity transported 
in the steady state (10). The exponents in i(//, T) are different from the exponents in 
L{p,Tqq) obtained when homology is applied to an entire star (Cox and Giuli 1968). 
We do not consider the dependence on surface photon temperature, TqQ, but on central 
temperature Tc- 

For A = 1, 1/ = 4.24 and the core opacity laws we consider, some numerical values are 
given in Table 1. (The table also contains the exponents for the Bohm-Vitense opacity. 
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with u = 4 for the ppl cycle.) The exponents are insensitive to the temperature exponent 
z/ = 4.24 in the luminosity generation law, but are sensitive to the opacity law. For the 
OPAL opacity, we obtain £ ~ m^'^^, in good agreement with the value L ~ M^'^ for the 
ZAMS Sun in Figure 3. 

We are interested in how the temperature varies as function of luminosity generation, 
opacity and mean molecular weight, for fixed luminosity Lq. 

Tc^i^iVeU)'/' . (15a) 

For the OPAL opacity function, we obtain the differential relation: 

dlnTc = (0.215)rfln;u - (0.133)rfln£o - (0.0344) In + (O.167)rflnL0 , (156) 

showing how the central temperature in any compact radiative core must change with input 
parameters. The central temperature is most sensitive to the chemically evolved mean 
molecular weight and to the overall luminosity generation Eo, and much less sensitive to the 
opacity, Kq- This is expected, since the core structure is determined by mass conservation, 
hydrostatic equilibrium and extended luminosity generation, while the radiative envelope 
structure is determined by the radiative transport and the central concentration of mass 
and luminosity. 

Because the energy generation is principally proportional to the ppl nuclear cross 

section factor, Eq ~ •S'n, we obtain Tq ~ S^^'^'^'^, in agreement with Iben (1969; 1991) and 

Castellani et al. (1993), who, however, incorrectly assumed p ~ T^. This explains why 
in Figures 2, 4, and 6-8 of Hata and Langacker (1995), Tc-parametrization is equivalent 
to parametrization, within the uncertainty of either. The ppl I, ppIII, and CNO 
chains contributions to luminosity generation break this simple Tc form, adding a weak pc 
dependence. The density exponents are: 

with ^ 

iP = a = ^ = 



o ^1 + ^ + ^ (17a) 
3 — v + s ^ ^ 

^ " 1 + X + n ' 

a = i{;- (er/e) 
b = a-iCT/9) 

d = T/e . 

The cases of the Kramers and OPAL opacities are presented in Table 2. The Bohm-Vitense 
opacity, with = 4, is also shown for comparison. 

A rotating core would change the hydrostatic equilibrium by adding the centrifugal 
force to that of gravity in the rest frame of solar matter. If, as in the magnetic case, 
the reverse influence of the thermomechanical structure on the rotation is neglected, the 
correction to the homology is: 

rVP~/m2[l-a;]3 

3(l + A)(3 + 2n) , 
^ 2(l + A + n) 
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where Iq is the non-rotating luminosity function, uj = Vp'r'^ /Gm, and O is the angular 

rotation frequency. Near the center, oj — > 30^/47rGpc- For typical SSMs, uJq ~ 2xl0~ , 
using a reasonable solar core rotation rate (Elsworth 1995). The exponent x is given 
in Table 2 for the three opacities. The rotation correction would be significant only for 
rotation rates ~ 400 times those in the Sun. Such a high rotation rate would not only 
change the thermomechanical structure, however, but could also induce sufficient chemical 
mixing to invalidate the homology. 

An alternative to this homology is to treat the luminosity not the radius r or the 
cumulative mass m, as the independent variable. Since the luminosity is a monotonically 
increasing function of r in the core, but not outside, this change of variables is feasible 
only in the core and separates out the luminosity-producing regions. 

4. Core Homology and Neutrino Fluxes 

After nuclear cross sections are introduced, and the ^He, Be abundances are assumed 
to be in steady state, each of the neutrino emissivities, fiy{i), is a function of X, p, T. 
Homology would then make each neutrino flux ~ fuii), subject to the luminosity 
constraint (1). If there were only one power-law energy generation term in equation (1), 
the core homology would be exact and pc, like all other core variables, would be a power 
of Tc only. The Be, B, and CNO neutrino production breakss this homology, so that, 
besides the principal sensitivity to Tc, the neutrino fluxes have a mild separate dependence 
on Pc- Using the luminosity constraint, Gough (1994) has obtained: 

HBe) ~ p^ ■ t! ~ p2-5^ ■ <PiBf'' (19) 
4>{B) ^ pP . ^ p-1-33 . 0(Se)2-33 . 

If we approximate p ~ in the solar core, we obtain ~ T^^^\ with a{pp) = 
— 1, a{Be) = 11, a{B) = 22; while CastcUani et al. (1993), assuming an n^fj = 3 poly- 
tropic Sun, obtained a(pp) = —1.1, a{Be) = 11, a(-B) = 27. The small departure from 
core homology, together with uncertainties in the nuclear cross-section factors 534, Sn, 
explains the scatter in diagrams plotting neutrino fluxes against Tc alone (Hata and Lan- 
gacker 1994; 1995). 

5. Core Homology and Helioseismology 

Helioseismology, the study of sunquakes, is based on three distinct types of waves in 
the solar medium, p-modes, /-modes, and g-modes (Hansen and Kawaler 1994). The flrst 
two are acoustic, with pressure contrast as the restoring force, and are seen in the outer, 
convective zone, where they have much or most of their amplitudes. Their eigenfrequencies 
rise with the number of nodes in the successive modes. 

The ^-mode restoring force is gravity, and these modes have their largest amplitude 
in the core. Their eigenfrequencies decrease with the number of nodes. The (7-modes have 
not yet been flrmly detected by optical means, although they have perhaps been detected 
through their modulation of the solar wind (Thomson et al. 1995). All modes are labelled 
by eigennumbers n and Z, with an azimuthal m if rotation is present. (Otherwise, the 
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eigenfrequencies are degenerate in m.) The ranges are: n = 1, 2, and / = 0, 1, 2, 
For large n, the frequencies of the ^f- modes are given by (Hansen and Kawaler 1994): 



f^^^^m ^ (20a) 
Jo ^ 

Nir) = ^9ir)/Xgir) , 

where g{r) is the local acceleration of gravity and Xg{r) a local scale height: 

1 1 11 1/ r\ 

. (206) 



^9 ^ad ^ Tad 

where Fg^d = ((ilnP/(ilnp)a^(j is the adiabatic polytropic exponent. 

Solar core homology can be applied to the integral Qg, which receives its main contri- 
bution from the core region. As r — > 0. 



^9^^d^(i;-^)^-Rc-mc)) , (21) 



where Rc is the core radius {Rc = {0.26)Rq) and 

is the radially averaged mean density interior to Rc- Note Fg^^ > F implies convective 
stability of the core. 

Because, outside the immediate central region {x < 0.049), the mass m(r) rises more 
slowly than r^, the integral emphasizes the central core as the dominant "yolk in the 
egg" mass concentration that controls the g-mode oscillations. A simple estimate is 
{p{Rc)) ~ Pc, but a better estimate results from applying core homology via the UqQ — oo 
"isothermal" polytropic solution (Chandrasekhar 1939; Kippenhahn and Weigert 1990). 

The important length scale here is ^fP^J^KGpl^ = O.O49i?0. Using the power series and 
asymptotic properties of the solution, one obtains: 

3pc f^" , dw{z)/dz 



{-p{Rc)) - ^ r dz 

Jo 

~ 0.56pc 



z (23) 



where Zc = 0.26/0.049 = 5.3 and w{z) is the dimensionless gravitational potential. This 
estimate is smaller and more accurate than pc as it covers the entire core, whose average 
density is lower than its central density. 

As helioseismological observations are so far in good agreement with SSM predictions, 
we conclude that the homology presented in this paper is, for practical purposes, a complete 
parametrization for any "reasonable" changes to the minimal SSM. 
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6. Conclusion 

Assuming mechanical and thermal stasis and neglecting chemical evolution, the homol- 
ogy makes the (Tc, pc) parametrization a general framework characterizing the inner core. 
Using this approach, we have disposed of three misconceptions: (1) that the luminosity- 
generating core of the Sun is polytropic; (2) that the polytropic relation p ~ T'^ is essential 
to understanding the Sun's core; and (3) that homology is inapplicable to stars on the 
middle of the main sequence. While the T-gradient depends on the opacity, Tc depends 
mainly on the mean molecular weight /i because of the homology, assuming a quiescent 
non-convective and unmixed core. The properties of the core depend only on one surface 
boundary condition, the total luminosity Lq, assumed to be in steady state with the core's 
luminosity. The surface temperature TqQ is then irrelevant. 

The solar core is almost homologous because its luminosity generation is dominated 
by the ppl cycle and, over its narrow range of temperature and density, the opacity and 
luminosity generation can be approximated by power laws. The luminosity of solar models 
based on purely radiative transfer scales by the nuclear cross-section factor Sn, or equiva- 

lently, by Tc ~ S^^'^"^. This quasi-homology justifies the Tc-parametrization as reasonable 
for estimating astrophysical uncertainties in any SSM and for extrapolating from any SSM 
to even extreme non-standard "cool Sun" or "hot Sun" models. In particular, analyses 
such as those of Hata and Langacker (1994; 1995) and of Bludman, Hata, Kennedy, and 
Langacker (1993), using the Tc-parametrization and nuclear cross section uncertainties 
alone, can arrive at theoretical neutrino fliix and detection rate uncertainties and their 
correlations, agreeing with those Bahcall and Ulrich obtained from 1000 different Monte 
Carlo SSM simulations (BahcaU and Ulrich 1988). 

Because of the central concentration of mass and luminosity generation, the Sun's 
outer radiative zone is nearly an n^Q = 4.3 partial polytrope, with exponential pressure, 
density, and temperature profiles. The density scale height Xp = O.O95i?0 is the single solar 
parameter entering into the MSW adiabaticity parameter that determines any small-angle 
MSW osciUations. 
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Appendix: The Luminosity Constraint 

The constraint of the total photon luminosity Lq upon the neutrino fiuxes is almost 
independent of the specific SSM, resulting for the most part from the microscopic properties 
of the nuclear fusion reactions (Schwarzschild 1958; Cox and Giuli 1968; Bahcall and Ulrich 
1988; Turck-Chieze 1988). The branching ratios assumed are mildly SSM-dependent and 
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are taken from the results of Bahcall and Pinsonneault (1992), with hehum diffusion. This 
luminosity constraint on the neutrino fluxes should be distinguished from that used to 
determine the helium abundance Y. 

The Sun shines by two nuclear reaction chains, the dominant pp and the minor CNO. 
The pp chain itself consists of three subchains, ppl, ppll, and ppIII, each terminating in 
'^He or a production in a different way. 

T 

pp/pep — > 

i 

I 
8S 



2a 



2a 



For each reaction, the energy released, Q, is partitioned between neutrino energy Qi^ and 
photon luminosity Q^, so that Q = + Qjy. Reactions without neutrinos have Qjy = 0. 

The luminosity constraint arises from the proportionality of the neutrinos fluxes to the 
nuclear reaction rates. Each reaction contributes its Qj value to the photon luminosity. 
Since there are two neutrinos emitted for each reaction chain, the luminosity is related to 
the neutrino fluxes by 

^0=(E<37W</'W)-(47ri?|/2) , (Al) 

i 

summed over all nuclear reactions i. Normalizing to the ppl value and to the neutrino 
fluxes measured at the Earth's orbit. 



Q^ippI) ^ttRI ' Q^ippI) ' ^ r, 



■0 

where = 149.6x10^ km is the Earth's average orbital radius. 

The Q, Qjy, and values for all four reaction chains are listed in Table 3 (Fowler 
1967; Bahcall and Ulrich 1988; Turck-Chieze et al. 1988). The pp chain is initiated by 
either the pp or pep reactions, in the ratio 99.6% : 0.4%; the ppll chain emits neutrinos 
at two discrete energies, in the ratio 89.7% : 10.3%. The quoted energies are weighted 
averages. The ratios of the values can then be computed to obtain the constraint (1). 

Also necessary are the ratios of the different neutrino fluxes in order to obtain equa- 
tion (2,3). These fluxes and the percentages above are obtained from a speciflc SSM, but 
the coefficients in the constraints (1-3) are relatively insensitive to theoretical variations 
that do not depart radically from conventional SSMs. The ratios of the reaction subchain 
termination rates are related to the flirxes by: 

term(pp//) 4>{Be) 



term(pp/) [(j){pp) /2] - (t){Be) ' 
teruijjrpIII) _ (j){B) ^ ' ' 

teTm{ppII) (l>{Be) 
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TABLE 1 



Luminosity Power Laws^ 


Opacity n 


s 


5 e C 




e 


Kramers 1 


3.5 


5.44 0.833 0.167 


1.33 


6.12 


BV 0.5 


2.5 


4.83 0.8 0.2 


1.3 


5.8 


OP AT n A'i 


2.47 


4.81 0.794 0.206 


1.29 


5.99 


^Exponents in 


the luminosity power-laws (15) for three Rosseland mean opacities of 


the form (7). 














TABLE 2 










Density Power Laws^ 






Opacity 


T 


a b c 


d 


X 


Kramers 0.333 


0.753 


0.231 0.313 0.497 


0.123 


2 


BV 0.4 


0.6 


0.317 0.372 0.534 


0.103 


1.8 


OPAL 0.412 


0.507 


0.345 0.395 0.521 


0.0847 


1.77 


^Exponents in 


the density power-laws (16) for three Rosseland mean 


opacities of the 


form (7). See text for last entry (x)- 










TABLE 3 










pp Reaction Q Values^ 






Chain 


Q 




Qj 




PPI{PP) 


26.732 


0.540 


26.192 




ppII{Be) 


26.682 


1.083 


25.599 




ppIII{B) 


26.639 


6.980 


19.659 





All energies in MeV. Qjy is spectrum-averaged for each reaction chain. 
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Figure Captions 



Fig. 1. The stiffness coefficient din P/ din p across the Sun's profile. In the convective zone, 
r = 5/3; in the outer radiative zone, F ^ 1.23; but over the core, F ~ 1 and varies. 

Fig. 2. The temperature gradient V = dlnT/dlnP and 1 — F~^ — dln{P/ p) /dlnP. Where 

the chemical composition is homogeneous, V = 1 — F"-*- = 0.4 in the convective zone 
and 0.18 in the outer radiative zone. But over the core, P/{lftpT) = (1 + D)/ p, is 
varying. 

Fig. 3. The lines shows calculated L-M and R-M relations for a large range of zero-age 
Main Sequence stars. The dots and triangles show best measurements of sleeted Main 
Sequence components of detached and visual binary components, respectively. (From 
Kippenhahn and Weigert (1990), by permission.) 
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r=dllnP/d!lnp 




